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ON AN EXPEEIMENTAL MODIFICATION OF VAN DEB 
WAALS’S EQUATION. 

By John P. Dalton, M.A., D.Sc., Victoria College, Stellenbosch . 


(Beceived October 11th. Bead October 15th.) 


§ 1. Introduction. 


Some years ago* the author published a complete table of 
saturation constants derived from the well-known van der Waals 
equation— 



(i) 


and, by their means, investigated the reliability of the equation in repre¬ 
senting experimental results, especially within the saturation region. 
The equation was found, as expected, to give an excellent qualitative 
indication of the nature of the results to be expected under given 
conditions, while it completely failed to give any reliable quantitative 
estimates. Many modifications of the above equation have been proposed 
from theoretical considerations based upon assumptions more or less 
convincing, but the arbitrary nature of such assumptions seems to 
suggest that a direct appeal to experiment is to be preferred. In the 
following an attempt is made on the lines indicated by Kuenen| to 
modify the equation experimentally, and the improvement effected is 
brought to light by comparison with such data as are at present 
accessible. 


* J. P. Dalton : Phil. Mag. (6), 13, 1907, pp. 517-542. 
f J. P. Kuenen : Die Zustandsgleichung (Vieweg), 1907, p. 128. 
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§ 2. The van der Waals “ a ” as a Function of the Temperature. 


Following Kuenen, we shall, as a first approximation, regard the a of 
the equation— 

(p+£)(*-*)-BT.(2) 

as a function of T only, and b as a quantity which is independent of T. 
By applying well-known thermodynamical relationships we easily 
obtain— 

. (3) 


where p s is the vapour-pressure at temperature T, and v T and v 2 are the 
specific volumes of the respective phases of the complex on the same 
isothermal. Integration between any two temperatures gives— 




W/ T 
T 2 V 



As experimental substance Isopentane was selected, for which the 
saturation data have been determined with great accuracy by Young,* 

"cm. 6 mm. Hg.~l 


P"" c 

From these data the integrand of (4) was evaluated in 


>C. 


units. These figures are given in column A of Table I. Integration was 
performed by means of a large-scale drawing on millimetre paper. The 
critical temperature was taken for the lower limit, and for the integration 
constant was written— 


fa \ 3 'p k v k 2 

T 


( 5 ) 


For isopentane this becomes— 

3 x 25015 x (4-266) 2 nnp , fern. 6 mm. Hg. "1 

-I«F8-= 2964 L'— * 0 — J 


Column B of Table I. gives the values obtained for the right-hand side 
of (4), and in the other columns are given values of a and of log a as 
functions of T. 

The last column of the table and Fig. 1 show conclusively that log a 
is a linear function of T, so that we may write— 


a = e a-/3T . 

* S. Young : Proc. Phys. Soc. 13 (1895), p. 602. 


(6) 
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Table I. 

a =f (T) : Isopentane. 


T°. 

A. 1 

1 

B. 

a 

T 

283 

48-37 I 

4493 

7457 

293 

44-59 i 

4066 

7030 

303 

41-07 , 

3674 

6638 

313 

38-00 1 

3310 

6274 

323 

35-30 I 

2970 

5934 

333 

33-11 

2653 

5617 

343 

30-76 

2357 

5321 

353 

28-92 

2079 

5043 

363 

26-87 

1821 

4785 

373 

25-11 

1578 

4542 

383 

23-62 

1351 

4315 

393 

22-09 

1137 

4101 

403 

20-74 

937 

3901 

413 

19-44 

749 

3713 

423 

18-26 

573 

3537 

433 

17-06 

409 

3373 

443 

15-76 

257 

3221 

449 

14-97 

— 

— 

453 

14-39 

106 

3070 

456 

13-83 

_ 

— 

458 

13-40 

— 

— 

459 

13-18 

'A- ' j 

— 

460 

12-89 

U. * * 

— 

460-4 

12-60 

1 — 

— 

460-8 

12-38 

0 

2964 


a x 10- 6 

log IC A 

log £ a. 

2-1103 

6-3244 

14-5625 

2-0598 

6-3138 

14-5381 

2-0113 

6-3035 

14-5143 

1-9637 

6-2931 

14-4904 

1-9167 

6-2826 

14-4662 

1-8705 

6-2720 

14-4418 

1-8251 

6-2613 

14-4172 

1-7802 

6-2505 

14-3923 

1-7369 

6-2398 

14-3677 

1-6942 

6-2290 

14-3428 

1-6527 

6-2182 

14-3179 

1-6117 

6-2073 

14-2928 

1-5721 

6-1965 

14-2680 

1-5335 

6-1857 

14-2431 

1-4962 

6-1750 

14-2185 

1-4605 

6-1645 

14-1943 

1-4269 

6-1544 

14-1710 

1-3907 

6-1432 

14-1452 

1-3658 

i z 

! 6-1354 

1 14-1273 
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The average value of /3 is 0-002446. This gives a = 15-255 in the units 
employed. 

With a view to the applicability of the law of corresponding states to 
the modified equation, (6) may be written in the form— 

. (7) 

where S is the reduced temperature, and /3T ft = lT27. 

§ 3. The Modified Equation. 

Equation (2) now takes the form— 

(p + ^ T * <I_ *')(»-&) = RT . (8) 

If the law of corresponding states is to hold good, /3T, C should be the same 
for all substances. Whether this is the case or not has not yet been 
determined, but it seems probable that it will not be found strictly 
constant, but to vary more or less with the complexity of the molecule, 
just as in the case of the reduced diameter slope.' 1 '- In these circumstances 
we may, as a first approximation to a general equation, take /3T /C to be 
unity, and equation (8) then yields on reduction— 

( 7r+ ^ £l "‘ & )( 3w ~ 1 ) = 8 $ . (9) 

In the sequel it is shown how much more closely this equation agrees 
with experiment than the original unmodified form. 


§ 4. The Modified Saturation Constants. 


A direct calculation of these constants is best made by a substitution 
method. If we write— 


3 Wl —1 = re~P ) 
3w 2 —1 = 7'e/ 3 i 


( 10 ) 


and substitute in (9) and in the Maxwell equation — 


irdit) = Tr s {ix) 2 — w r ) 


we can eliminate 7r and S and obtain- 


(ii) 


n £ 2 P—£~ 2 P — 4/3 

^ = (/ 3 + l)£-P+((i-l)£fi 
sinh 2/3 —2/3 
~ j3 cosh /3 — sinh (3 


( 12 ) 

(13) 


* Cf. E. Mathias, H. Kamerlingh Onnes, and 0. A, Crommelin : Comm. Leiden, 
No. 181a (1918). 
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r is now evaluated from tables for various values of the parameter /3, and 
the saturation constants then obtained. Even with tables the calculations 
are somewhat laborious, and it was found to be much easier to derive the 
modified constants from the table of unmodified constants referred to 
above. It is evident that the modification now introducted does not 
affect in any way the relative values of the saturation volumes ; * hence 
if TI and 0 are the ordinary van dee Waals saturation pressure and 
temperature for any pair of volumes, and 7 r and $ represent the modified 
constants, we find— 


and— 


1 e 1 -^ 

e = T~ 


(14) 


n 7r 

Q = f ' 


(15) 


If equation (14) is written in the form — 


log c 100 = log e 10 $ — 1 + $ 

it can be solved by inspection from a table of natural logarithms ; 7 r is 

then obtained from the known ratio Proceeding in this way, the 

0 

following table of modified constants was obtained. 


Table II. 

Saturation Constants according to the Modified Equation. 


S’. 

7T. 

| ... 

w 2 . 

3 t . 

S 2 . 

0-4964 

0-0005275 

0-3698 

2506 

2-704 

0-0004 

0-5493 

0-002460 

0*3777 

592-6 

2-646 

00017 

0-5980 

0-007736 

0-3864 

203-6 

2-587 

0-0049 

0-6430 

0-01877 

0-3960 

89-15 

2-525 

0-01122 

0-6850 

0-03807 

0-4068 

45-98 

2-459 

0-02175 

0-7245 

0-06795 

0-4188 

26-61 

2-388 

0-03758 

0-7615 

0-1103 

0-4326 

16-73 

2-312 

0-05977 

0-7966 

01665 

0-4485 

11-18 

2-230 

0-08946 

0-8298 

0-2376 

• 0-4672 

7-811 

2-140 

0-1280 

0-8614 

0-3244 

0-4896 

5-643 

2-042 

0-1772 

0-8916 

0-4272 

0-5174 

4-173 

1-933 

0-2396 

0-9204 

0-5463 

0-5534 

3-128 

1-807 

0-3197 

0-9480 

0-6815 

0-6034 

2-349 

1-657 

0-4257 

0-9745 

0-8328 

0-6841 

1-727 

1-462 

0-5790 

1-0000 < 

1-0000 

1-0000 

1-000 

1-000 

1-000 


If, as is in the meantime assumed, 6 is a constant. 
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The improvement effected by this modification is strikingly shown in 
Fig. 2. The reduced vapour-pressures from equations (1) and (9) are 
plotted, together with the experimental curves for isopentane * and methyl 
alcohol. If the law of corresponding states were rigorously obeyed the 
last two curves would be coincident; as long as such differences as these 
occur it is obviously impossible to represent the thermal behaviour of 
different substances by a single equation, and all that one can do is to 
construct an equation for a fictitious substance to serve as a type from 
which particular substances deviate more or less according to their 
molecular complexity. This is what Kamerlingh Onnes has done in 
deducing his “mean empirical reduced equation.” Equation (9) seems, 
at least in so far as vapour-pressures are concerned, to be as good a type 
as can be found to represent both normal (isopentane) and abnormal 
(methyl alcohol) substances; it gives a much closer quantitative estimate 
than the unmodified equation, while, at the same time, it retains all the 
essential features of the original. 

From the following table it will be seen that the modified equation is 
in good agreement with the van der Waals empirical vapour-pressure 
law— 


log*=/(l-i).(16) 

as long as the temperature is not too low, and that the value of / at the 
critical point, 3'04, is not far from the corresponding values for such 
typically normal substances as carbon dioxide, 2-97, and isopentane, 2-95. 
The unmodified equation, on the other hand, as was previously shown,! 
gives values for / which deviate greatly from constancy, while the critical 
value is 1-72. 


Table III. 


3. 

/• 

3. 

/• 

3. 

/■ 

0-496 

3-23 ; 

0-725 

3-07 

0-892 

3-04 

0-549 

3-18 

0-762 

3-06 

0-920 

3 04 

0-598 

3-14 , 

0*797 

3-05 

0-948 

3-04 

0-643 

3-11 1 

0-830 

304 

0975 

3-04 

0-685 

3-09 

0-861 

3-04 

1-000 

— 


* The carbon dioxide carve is coincident with this, 
f Loc. cit. 
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§ 5. Latent Heat. 

As a further test of the improvement effected by the foregoing 
modification, we may use the new constants for calculating values of 
the latent heat at different temperatures. The external work (reduced) of 
vaporization— 

W .(17) 

is only a fraction of the total latent heat, and yet even there an improve¬ 
ment is already evident, quite marked with respect to dependence on 
temperature although the change effected in the absolute magnitude of the 
work is not so great. Van der Waals’s equation gives 1'42 at 3 = 0'6, 
rising to a maximum of L47 at 3 = 0‘7 and, of course, falling to zero at the 
critical point; while at 3 = 061 isopentane gives 2 - 22, rising to a maximum 
of 2 - 42 at 3 = 076. Values calculated from the saturation constants of 
Table II. give 157 at 3 = 060 rising to a maximum, L80, at 3 = 0-76. As 
before, carbon dioxide * and isopentane correspond very well. 

The experimental values of the latent heat have been obtained from 
the well-know T n equation— 

L = T^a-*,).(18) 

and the following table contains the values of the reduced latent heat, 
L>/p k v k , thus found for isopentane. 

Table IV. 


-- : Isopentane. 


3. 

L 

PkVk 

3. 

L 

PkVk 

3. 

L 

PkVk 

0-614 

25-00 

0-788 

19-88 

0-961 

10-09 

0-636 

24-41 

0-809 

1907 

0-974 

8-60 

0-658 

23-72 

0-831 

18-19 

0-983 

7-257 

0-679 

23-17 

0-853 

17-23 

0-990 

5-991 

0-701 

22-59 

, 0-875 

16-21 

0-994 

4-849 

0-723 

21-95 

0-896 

1503 

0-996 

4-127 

0-744 

21-22 j 

0-918 

13-74 

0-998 

3-071 

0-766 

i 20-61 

; 0-940 

12-15 

0-999 

2-376 


* J. P. Kuenen and W. Gr. Robson ; Phil. Mag. (6), 3, 155 (1902). 
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If van der Waals’s equation is used for the calculation in con¬ 
junction with the thermodynamical relation— 

/*V 2 

L = tJ .(19) 

we find— 


L 8 , 

^“3 el °g- 


3w 2 —1 
3w x — 1 ■ 


( 20 ) 


while the modified equation gives— 


L 

PkV k 


—(H) • ■ <*> 


When numerical values are obtained for (20), by means of the saturation 
constants of equation (1), we get, as the following table shows, reduced 
latent heats which are exactly one-third of the isopentane values at 
the same temperature. The same holds for carbon dioxide, for the latent 
heat curves for C0 2 * and for isopentane practically coincide when 
reduced. 


Table Y. 


-: VAN DER WAALS. 

p k v k 


3. 

L 

PkV k 

3. 

L 

PlPk 

1 

3. 

L 

PkVk 

030 

8-91 

0-55 

8-40 

0-80 

6-48 

0-35 

8-87 

0-60 

8-17 

0-85 

5-76 

040 

8-80 

065 

7-88 

0-90 

4-82 

0-45 

8-71 

0-70 

7-51 

095 

3-50 

0-50 

8-58 

0-75 

7-05 

1-00 



When the latent heat is calculated from (21) instead of (20), the 
greater part of this large discrepancy between actual and calculated 
values disappears. In the table, column A gives values of the first 
term of (21) and B gives values of the second term, each being obtained 
from the saturation constants of Table II. 


Kuenen and Robson : loc. cit. 
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Table VI. 


— - : Modified Equation. 



A. 

B. 

L 

PkVk 

0-496 

14-75 

6-61 

21-36 

0-549 

13-92 

6-84 

20-76 

0-598 

13-16 

6-92 

20-08 

0-643 

12-44 

6-90 

19-34 

0-685 

11-75 

6-86 

18-61 

0-725 

11-07 

6-73 

17-80 

0-762 

10-37 

6-53 

16-90 

0-797 

9-66 

6-28 

i 15-94 

0-830 

8-91 

5-94 

14-85 

0-861 

8-10 

5-54 

13-64 

0-892 

7-22 

5-05 

12-27 

0-920 

6-24 

4-45 

10-69 

0-948 

5-08 

3-69 

8-77 

0-975 

359 

2-64 

6-23 

1-000 

0-00 

0-00 

0-00 


§ 6. On the Specific Heat of Saturated Vapours. 

In one of the papers to which reference has already been made (§ 1) 
the present writer used van der Waals’s equation to investigate the 
behaviour of the specific heat of saturated vapours as a function of the 
limiting value of the ratio of their specific heats. As data were then (and 
still are) lacking by means of which the conclusions reached could be 
tested directly, an endeavour was made to obtain an experimental curve 
from the isopentane isothermals under a certain assumption regarding the 
behaviour of the specific heat at constant volume. That assumption has 
been criticized * as being not in accordance with the law of corresponding 
states, and it has, therefore, been abandoned in the following treatment 
of the problem. Specific heats are first obtained by means of the modified 
equation (which has been shown to approximate more closely to ex¬ 
periment than the van der Waals equation), and an experimental 
curve is then obtained from the isopentane isothermals, using the 
Kamerlingh Onnes equation of state for determining the change of 
with volume. 

* H. Kamerlingh Onnes and W. H. Keesom : Die Zustandsgleichung, Encyc. Math , 
Wiss., v. 10, 987 (1912) or Comm. Leiden. Suppl. No. 28, p. 323, 
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If G s represents the specific heat of the saturated vapour we may write 


r _p , t _ 5 p *>• 

o.-o,+ i aT dT 


while C, is connected with C Ur by the relation — 


( 22 ) 


C„ = C„,+T| 


(23) 


We may also write— 

. w 

Combining these results we find for the specific heat of the saturated 
vapour in terms of reduced magnitudes — 


C f - 


\ rh, i ^dta, ) / 9 k\ 

1= tA i av rf " + s1 ■ ■ (25) 


730 - 

Along the inversion curve C s vanishes; its equation is therefore given by— 

8 


= 1 ■ 


From equation (9) we find- 


q ft | f d a 7T J , ^7T dh) s j 

35 ) + } 
Jo 


• . (26) 


07T _ 8 

~ 3oJ=T ' 


3 J-S 


■ (27) 


and— 



so that the inversion curve becomes — 


(28) 


7 


co 


5? e l —7 _i_ d^. s f 3 | 

8 oi ^ 



(29) 


The right-hand member of the equation was evaluated from the constants 
of Table II., and the following results were thus obtained :— 
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Table VII. 

Inversion Curve: Modified Equation. 


&. 

7oc • 

3. 

7oo • 

3. 

7 00 • 

0-496 

1-066 

0-725 

1-099 

0-892 

1-083 

0-549 

1-075 

0-762 

1-100 

1-920 

1-073 

0-598 

1-083 

0-797 

1-099 

0-948 

1-060 

0-643 

1-090 

0-830 

1-096 

0-975 

1-043 

0-685 

1-095 

0-861 

1-090 

1-000 

1-000 


To obtain the experimental inversion curve, isopentane was again taken 
as a typical normal substance and equation (26) was applied to it in the 
unreduced form— 


• 1 -- 


*p dv , T i2_*i£ 


(26') 


The second term of the denominator has already been evaluated by 
Dieterici* for isopentane, and his values have been used. The other 
term, which occurs only as a small correction, was evaluated by means of 
the Kamerlingh Onnes mean reduced equation.! According to that 
equation the product, 7 rw, is expressed as a finite series of inverse powers 
of the reduced volume— 


,ry = $(i + B- + C-+D*± + ...) .... (30) 

K \ (O o ) 2 w 4 J 

where B, C, D, etc., are expressed as polynomial functions of the 
temperature. The complete equation has five such “ virial-coefficients,” 
but for the present purpose the three given above were sufficient. 
BT 

k is the critical constant - - 

P*v» 

We now find, to the degree of approximation required— 




- ll(BS)+- ll(C.9)+ -Jt (DJ) 
io 2 dy K 1 W 3 ds* y ’ oj5 dS 2 v 


(31) 


and from this we easily obtain the correctional term— 



<**>+£ s 33 =< OS)+ 


fa d 2 

4 E~^dS 2 



* C. Dieterici: Ann. d. Phys. 12, 602 (1908). 
f Kamerlingh Onnes and Keesom: op . cit . p. 729. 


( 32 ) 
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From the constants given by Kamerlingh Onnes and Keesom ( loc. 
cit.) this expression was evaluated along the isopentane saturation line. As 
the Dieterici values are a little irregular at the critical part of the curve, 
they were smoothed graphically before being used. In each case figures 
are given in cal/° C. and are for the gram-molecule. 


Table VIII. 

Experimental Inversion Curve. 


Sr. 


( smoothed )- 


0-636 

0-5 

-24-0 

108s 

0-679 

0-7 

-22-7 

1-09° 

0-723 

1-0 

-21-3 

1-09 8 

0-766 

1-3 

-19-8 

1-107 

0-809 

1-6 

-19-1 

l-lla 

0-831 

1-8 

-19-4 

1-113 

0-853 

1-9 

-19-8 

1-11* 

0-875 

2-1 

-20-6 

1-107 

0-896 

2-3 

-21-9 

1-10 1 

0918 

2-6 

-240 

1-093 

0-940 

2-9 

-28-2 

1-07 8 

0-961 

3-3 

-34-8 

1-063 

0-983 

3-8 

-53-8 

1-04° 

0-994 

4-4 

-114-4 

1-01 8 

1-000 

— 

— 

— 
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Both curves (Big. 3) exhibit the same general features as those 
previously obtained—for y^ greater than a certain critical value, C s is 
negative for all values of 9 ; for smaller values of y^, we have C s becoming 
negative, zero, positive, zero, and again negative as 5 increases towards 
unity; the two inversion temperatures are further apart the lower the 
value of y^ , and they coalesce at the critical y^ . But this critical value 
is somewhat smaller than expected, and, if the values previously quoted 
for y-Q are correct, would seem to exclude the possibility of inversion 
occurring with Acetone, Benzene, and Chloroform, for which it has actually 
been observed. But it is quite possible that the values of y which were 
then used were obtained under conditions not sufficiently far removed 
from the state of saturation to allow of their being regarded as identical 
with the state of a perfect gas. The Wiedemann y for Benzene, for 
instance, was 1-129, while a more recent determination by Stevens * 
gives (at + 100° C.) 1T05; according to the present curves the former 
value would render C s always negative, while the latter would give in¬ 
version at -9 = 0-75 and 9-= 0-88. For isopentane, for which y^ is 
probably 1-08, the experimental curve would give inversion at 9 = 0’61 
and 9 = 0-74, values which are, as one would expect, almost identical 
with those deduced by Dieterici from his own experiments in con¬ 
junction with Young’s isothermals. But the paucity of experimental 
data for both y^ and G s renders it impossible to test the quantitative 
value of the foregoing conclusions. 

October , 1913. 


E. H. Stevens: Ann. cl. Phys. (4), 7, 285, 1902. 


